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The nature of dark matter is one of the most thrilling riddles for both cosmology and particle 
physics nowadays. While in the typical models the dark sector is composed only by weakly in- 
teracting massive particles, an arguably more natural scenario would include a whole set of gauge 
interactions which are invisible for the standard model but that are in contact with the dark matter. 
We present a method to constrain the number of massless gauge bosons and other relativistic parti- 
cles that might be present in the dark sector using current and future cosmic microwave background 
data, and provide upper bounds on the size of the dark sector. We use the fact that the dark matter 
abundance depends on the strength of the interactions with both sectors, which allows one to relate 
the freeze-out temperature of the dark matter with the temperature of this cosmic background of 
dark gauge bosons. This relation can then be used to calculate how sizable is the impact of the 
relativistic dark sector in the number of degrees of freedom of the early Universe, providing an inter- 
esting and testable connection between cosmological data and direct/indirect detection experiments. 
The recent Planck data, in combination with other cosmic microwave background experiments and 
baryonic acoustic oscillations data, constrains the number of relativistic dark gauge bosons, when 
the freeze-out temperature of the dark matter is larger than the top mass, to be A*' < 14 for the sim- 
plest scenarios, while those limits are slightly relaxed for the combination with the Hubble constant 
measurements to N < 20. Future releases of Planck data are expected to reduce the uncertainty by 
approximately a factor 3, what will reduce significantly the parameter space of allowed models. 

PACS numbers: 98.80.-k, 98.70.Vc, 95.35.+d, 98.80.Cq 



I. INTRODUCTION 

The standard model (SM) of particle physics has been 
extremely successful in describing accelerator and other 
terrestrial experiments over the last decades. The recent 
discovery of a Higgs-like particle at the Large Hadron 
Colhder [HE], for instance, seems to piece together a 
model that describes the fundamental interactions of 
nature using gauge theories. Nonetheless, cosmological 
data indicates that the matter content described by the 
SM are responsible for only around 4% of the current 
energy density of the Universe the remaining be- 
ing in the form of the still unknown dark energy that 
drives the acceleration of the cosmos [4 , and dark mat- 
ter (DM) which holds together galaxies and larger struc- 
tures and corresponds to about 26% of the energy budget 
of the Universe. Dark matter is assumed to be com- 
posed by particles, and typical candidates are in the 
form of weakly-interacting massive particles (WIMPs) 
that arise naturally in supersymmetric models [5] , extra- 
dimensions [7] , and extended Higgs sectors E] . In 
analogy to what we observe in the visible sector of the 
Universe, the dark sector (DS) could be composed by 
more species, besides the DM candidate, including extra 
dark particles and interactions, as for instance the case 
of models with hidden SU(A^) gauge symmetries. This 
scenario would certainly lead to a richer phenomenology, 
if not in accelerators, at least in terms of the early and 
late Universe. Those species could, for instance, freeze- 
out/decouple at different temperatures and leave an im- 
print in the early Universe, and/or could potentially ex- 
plain the DM abundance in a different fashion than the 
so-called "WIMP-miracle" [lOl [II] where the DM abun- 
dance comes from the SM gauge and/or Higgs bosons 
exchange between DM and SM particles. Other alter- 



natives include the so-called mirror DM models [T^ [T5] 
where the interaction between DM and SM particles hap- 
pens through the kinetic mixing between SM photons 
and dark photons, the latter coming from a hidden U(l) 
gauge symmetry interacting only with DM. The strength 
of the mixing controls then the DM abundance, indirect, 
and direct detection signals. Moreover, this mixing al- 
lows dark photons to slightly interact with SM charged 
particles, e.g. with electrons and positrons, modifying 
the evolution of the Universe [14j. Nevertheless, such ki- 
netic mixing is not present when the DS is composed by 
unbroken non-abelian gauge symmetries. 

In this article we consider the implications that in- 
teracting WIMP-DS scenarios would have for the early 
Universe cosmology. We discuss how current cosmic 
microwave background (CMB) observations from the 
Planck satellite [H], in combination with the Atacama 
Cosmology Telescope (ACT) [l6], the South Pole Tele- 
scope (SPT) [H], WMAP data [H], baryonic acoustic 
oscillations (BAO) [T^22j . as well as with astrophysical 
measurements of the Hubble constant [231 El] , can place 
limits on those new light particles within the dark sec- 
tor. As we will show, if those "dark radiation particles" 
responsible for the DS interaction were in thermal equi- 
librium with the visible Universe via DM particles, they 
could potentially leave an imprint in the number of de- 
grees of freedom (d.o.f.) of the early Universe that could 
be observed and/or constrained in the near future. 

While similar scenarios have been discussed before in 
the literature (see for instance [25H28] and references 
therein), our discussion generalizes the calculations for 
potentially more realistic situations in which the interac- 
tions in the DS are mediated by more than one species 
that can have different freeze-out temperatures. More- 
over, our results do not exclude other cosmological effects 
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similar to the case of neutrinos in the visible sector, as 
the currently cosmological model requires at least three 
neutrinos with a temperature slightly lower when com- 
pared to photons, ~ (4/11)1/37:^ pg] . 

If those interactions take place within the dark sector 
only, cosmological observations like the ones discussed 
here will potentially be one of the few windows available 
to understand those new fundamental sectors. 



II. EARLY UNIVERSE COSMOLOGY 



Ti = XiTi) and masses uii <^ XiTo- The latter is in- 
spired by the visible sector, in which neutrinos and pho- 
tons have different temperatures. Tq is the temperature 
of the background of dark gauge bosons, which is lower 
than the temperature of visible photons as the latter get 
entropy contributions from particle-antiparticle annihila- 
tions. To calculate this extra contribution to N^Bi one 
can use the entropy conservation to estimate the ratio of 
temperatures in both sectors. Since the entropy is given 
by S' = g^si^T^, being a the scale factor, and under the 
approximation that all the decouplings occur instanta- 
neously, the ratio of temperatures can be written as 



The evolution of the Universe during the radiation 
dominated epoch is described by the standard Friedmann 
equation [2^1^. 
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where H is the Hubble rate and G the Newton constant. 
The energy density in radiation after the annihilation of 
electron-positron pairs is given by 



Pr = P7 



7 /II 
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X Nr, 
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Here the effective number of neutrinos is given by 
iVeff = + AA^cff , where N^f = 3.046 is the stan- 
dard contribution from the three SM neutrinos [?TH55] . 
and AA'cff is the contribution coming from extra light 
species present in the early Universe, known in the 
literature as dark radiation^ like e.g. sterile neutrinos 
or from other physical processes, like lepton 
asymmetries [351 EOj: reheating of neutrinos [3T], grav- 
itational waves [42], and other effects [43]. Here, for 
simplicity of the analysis, we will assume that those 
other possible contributions are negligible, and any extra 
contribution for TVeff comes from the extra d.o.f. stored 
in the dark sector. 
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where gl:°' is the number of d.o.f. of the Universe when 
the interactions that maintained the dark (relativistic) 
and visible sectors in thermal equilibrium ceased, and 
5*s0 corresponds to its current value. For thermal DM 
candidates, can be related to the DM freeze-out tem- 
perature and to the dark matter mass Mdm- In the 
case of WIMPs, for instance, we have that 
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M, 



DM 
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is a fraction of the DM mass [44] , and therefore one could 
go a step further and relate Td to Mdm, a feature that 
holds for any model in which the DM is the only link 
between the dark and visible sectors. In addition to that, 
the observed DM abundance [15 , 



r^D^/i =0.1196 ± 0.0031, 
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requires 
of the 



a precise value (to the percent level) 
thermally averaged scattering cross-section 



iav) 



3 X 10 cm^s ^, what fixes the relation among 



DM and its interactions with the visible and dark sectors. 



{(jv) = (cru)sM + (o-w)ds 3 X 10 



cm^s 



(6) 



III. DEGREES OF FREEDOM OF THE DARK 
SECTOR 

To account for the d.o.f. in the dark sector and cal- 
culate the contribution to the effective number of neu- 
trinos, we will adopt the standard tools for calculating 
the entropy and gravitational impact of the new parti- 
cles in the early Universe [3D]. We assume that those 
new particles are either massless or very light, and there- 
fore relativistic, when their interaction rate F^i in the 
dark sector became smaller than the Hubble rate 
effectively freezing-out their distribution as bosons or 
fermions. In our analysis the dark sector can be com- 
posed by N species of dark gauge bosons (DBs), i.e. 
massless spin-1 particles, that mediate the interactions 
between dark matter particles and M species of rela- 
tivistic dark fermions (DFs) with different temperatures 



and this, in turn, can be translated into a relation be- 
tween the dark sector decoupling and dark matter freeze- 
out temperatures [50] . 



7^f .0. 

-'dm 



18 + log 



{(jv) 



Ah 



DM 



3 X 10-27 cm^s-i GeV 
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3 X 10-27 cni3s-i GeV 
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The case when dark radiation particles cannot commu- 
nicate directly with the visible sector allows DM to help 
dark radiation particles to be in thermal equilibrium with 
photons up to the time of DM freeze-out. This condi- 
tion connects the temperature of DM freeze-out with 
and roughly implies a lower bound on Td set by 
when the DM annihilation cross-section is dominated 
by DS interactions, as can be seen from Eq. ([t]). This 
case, ((Tf)Ds — ('7''^) 7 would reduce the posibility of de- 
tection of WIMPs in direct/indirect searches. On the 
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other hand, a scenario with WIMP-detection prospects 
would be: (cr?;)DS ~ 10~^ x (cru), which increases the 
decoupling temperature to Td ^ 1.5 x due to the 
logarithmic dependence shown in Eq. ([7|. This scenario 
can be well motivated in the context of particle physics 
models where there exist a hidden gauge symmetry like 
SU(A'^) symmetries. These hidden symmetries can reduce 
the interaction strength between DM and the SM and, 
as consequence, the indirect and direct detection signal. 
This feature is important to avoid exclusion regions set 
by direct detection experiments like XENONIOO 45,. In 
a similar fashion, asymmetric DM models (that aim to 
explain why DM and baryon abundances are similar) can 
affect the evolution of the early Universe when interac- 
tions with the dark radiation particles are included |27j . 
However, contrary to the WIMP case, for those scenarios 
the mechanism of DM production and freeze-out is usu- 
ally model dependent, and therefore a general study is 
more challenging and probably not feasible. In the sce- 
narios discussed here, the main prediction is the connec- 
tion of the freeze-out temperature with the extra d.o.f. 
that will impact -/Voff as, for instance in eq. ([7|, that we 
will discuss in detail elsewhere 

Using the definition of entropy d.o.f. and eq.([3| one 
can write g■^,so as 



g*so ~ 3.91 + gs,D [ -j;^ 



where 



N M 
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(9) 



and therefore, 



9*s0 ~ 3.91 , f r 



(10) 



The increment in N^s due to the new particles in the 
dark sector can then be written as 



where 



3 91 

ATVeff « 2.201 <?D ( 

\.9*s - .9s,D 



N ^ M 

9B = Y^g-,xf + -Y^gjx'^ 



4/3 



(11) 



(12) 



Notice that the factor 3.91 comes from the number of 
entropy d.o.f. after the e'^e'^ annihilations. 



kind of constraints can currently be placed in those sce- 
narios using CMB data (in combination with other cos- 
mological observations). First, we will describe a sce- 
nario in which all TV massless "dark gauge bosons" have 
the same temperature, and it is associated with the dark 
matter mass via the freeze-out temperature of DM. For 
simplicity, we assume that all DBs belong to the same 
family of the dark particle background i.e. xdb.i = 1. Af- 
ter that, we will generalize this scenario by allowing for 
M massless "dark fermions" that have a slightly lower 
temperature than the dark bosonic sector, a situation 
that mimics the photon and neutrino backgrounds of the 
visible sector. 

In the case of N dark gauge bosons, eqs. ([9|, ([IT]), and 



( 12 1 connect the three main quantities as 



AiV, 



off 



4.402 iV 



3.91 
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In this case, glf' and N have a non-null impact on 
iVefr except in trivial cases like iV = or extremely large 
Ti i.e. glf' — > oo. The implications for N^g are presented 
in Fig. [l] where the lines show two of the 95% C.L. con- 
straints quoted in Ref. |T5l: from CMB (where CMB 
corresponds to the combination of data from Planck, 
WMAP-Polarization, SPT, and ACT) on N^g [ig, in 
combination with (i) BAG 19-22J, 



o on+0.54 



(95%C.L.;CMB-|-BAO), (14) 



and in combination with (ii) the measurement of the 
Hubble constant IHIll- 



iVoff = 3.62l^ j° (95% C.L.; CMB+i?o 



(15) 



Notice that the latter excludes at 95% C.L. the standard 
value of Ncs, requiring the presence of extra radiation, 
such as the dark sector. The diagonal (light dotted) lines 
represent the number of total degrees of freedom (visible 
and in the dark sector) at the freeze-out epoch, glf' = 

.o. I I .o. I 

9*s \sM + g*s Ids- 

Depending on the decoupling temperature of DBs, we 
can infer the maximum size of the dark sector allowed 
by the uncertainties in A'eff. For instance, all models 
for which is larger than the top mass "live" in the 
line indicated by T^^^ > rrit, and therefore current con- 
straints (dashed line) on N^e present an upper limit on 
the number of dark gauge bosons of approximately 14 
(CMB-I-BAO datasets). For temperatures between the 
bottom-quark and the W, Z bosons masses, the bound is 
reduced to « 10 DBs. Smaller values of shrink the 
dark sector down to few or no DBs, although we should 
keep in mind that sub-GeV temperatures are affected by 
the QCD phase transition. 

The addition of M = 3A^ dark fermions affects the 
value of Ncs by 



IV. DISCUSSION 
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We will discuss here two particular toy models that 
simplify the equations above with the goal to learn what 



which is similar in form to the case composed only by 
DBs. In Fig. [2] we present the implications for this setup. 



4 




5 10 15 20 25 

N {- number of gauge bosons) 



FIG. 1: Number of d.o.f. glf' versus number of dark gauge 
bosons. The red dashed line corresponds to the current 95% 
C.L. upper value on N^s using CMB+BAO [15], while the 
upper/lower blue solid lines shows the combination CMB+Ho 
(see text for details). Diagonal dotted lines represent glf' 
levels associated to Td for the SM particle content. Notice 
that for a given range of T^m, models exist only along the 
corresponding line. 



As one would expect, the presence of DFs leads to a re- 
duction on the maximum number of DBs because they 
also contribute to the radiation density. For the freeze- 
out temperatures discussed in the previous case, we ob- 
tain an upper bound of approximately 10 and 9 dark 
gauge bosons for > rrit and nib < ?dm < ^w,z, 

respectively for current data. The expected constraints 
that one can place on this scenario are expect to improve 
in the near future as CMB experiments continue to take 
data, and we expect that future data releases might rule 
out whole classes of models. 

We should emphasize that the dark sector also con- 
tains one or more stable (or quasi-stable) massive states 
that play the role of DM. While those particles have a 
negligible contribution to our calculations in terms of ra- 
diation, they have governed the evolution of the Uni- 
verse for a long time and are responsible for the forma- 
tion of structures. For simplicity, we consider the case of 
WIMPs, which is a widely studied DM candidate with 
well-behaved thermal history. If the DS interaction rate 
is similar to the SM interaction rate, Eq. ^ implies a 
relation between and Mdm- In particular, WIMP 
masses in the range 100 GeV - 1800 GeV constrain the 
maximum number of DBs to « 14 (20) at present, de- 
pending on the datasets. 

While this analysis was done under the assumption 
that the dark and visible sectors shared the same 
temperature at some point, it should also work for 
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FIG. 2: Number of d.o.f. glf' versus number TV of dark 
gauge bosons for the case with additional 3A^ dark fermions. 
Line-styles are equivalent to the ones shown in Fig. [l] 

models that were never in thermal equilibrium with SM 
particles but for which the couplings of the inflaton with 
both the dark and visible sector are known: in this case 
one can show that Td will be related to the reheating 
temperature and to the couplings of the inflaton with 
each species [25], and they can be translated into 
different Xi. Yet another alternative is to think of the 
temperature ratios Xi as free parameters instead of being 
given by the particle physics model, and use future 
precise data to constrain them, although in this case one 
would most probably need extra datasets to disentangle 
the correlations between Xi and the number of particles 
in the dark sector. 



V. CONCLUSIONS 

The dark sector is generally assumed to be composed 
only by a WIMP dark matter, but a feasible (and ar- 
guably more natural) model would include dark gauge 
bosons and other particles that might be stable. In this 
article we discussed what is potentially one of the very 
few windows to place constraints on the number and in- 
teractions of those particles, and applied our general re- 
sults to a couple of toy models to understand to what 
level current and future cosniological data can constrain 
those scenarios. 

In particular, we showed that for models in which 
the DM is weakly interacting, one can relate its freeze- 
out temperature to the temperature of the rest of the 
dark sector, and therefore amend the standard analysis 
of d.o.f. in the early Universe to take into account the 
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contribution of this dark sector to the effective number 
of neutrinos. Moreover, in the context of the standard 
cosmological model, our method can also be easily gener- 
alized for cases in which the dark sector was not in ther- 
mal equilibrium with SM particles, since both sectors are 
supposed to have a common origin after the inflationary 
phase. 

We calculated this extra contribution to N^s in an al- 
most model independent way, and particularized those 
results for some simple models to compare them with 
current and future data. In particular, we obtained the 
constraints on the number of interactions within the dark 
sector: while current limits on the number of dark gauge 
bosons are somewhat weak (< 14 (20)), future CMB data 
will potentially improve those bounds by a factor of ap- 
proximately 3 (assuming they can improve sensitivity to 
f [iVoff] « 0.2 at 95% C.L.), ruling out whole classes of 



models. It will be interesting to look at other WIMP DM 
models that can be accommodated within this scenario 
and constrain their properties with current and forth- 
coming cosmological data. 
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